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Abstract 

We report on a microscopic Refined Resonating Group Model (RRGM) calculation 
of scattering of p off 3 He employing the Argonne--ui4 and the Bonn nucleon-nucleon 
potentials without three-nucleon forces at low energies up to 30 MeV. The calculated 
phase-shifts verify the well-known proton analyzing power A^-problem. We demon- 
strate that with corrected 3 P2 phase-shifts experimental differential cross-section 
and analyzing power data can be explained. 
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1 Introduction 



It is widely known that realistic nucleon-nucleon (NN) forces cannot repro- 
duce the 3 H and 3 He binding energies. The A y analyzing power in (N — d) 
scattering is still not reproduced by adding three-nucleon interactions [1]. The 
30% deviation of A y can be resolved by tiny changes on the order of 0.1° 
in the (N — d) scattering phase-shifts [2,3,4]. It was shown in [5] that al- 
though a realistic NN-force can generally reproduce the 4 He system, there 
remain differences, most notably in the analyzing powers. Recently Fonseca 
et al. [6] studied the (n — 3 He) scattering system using an AGS-approach. 
They calculated the cross-sections and vector analyzing powers using the NN- 
interactions Bonn-CD [7], Nijmegen II [8] and Bonn-B [9] without 3A-forces 

Email address: reissc@kph.uni-mainz.de (C. Reifi). 



Preprint submitted to Elsevier Science 



8 February 2008 



and pointed out that the (n — 3 He) system may turn into a A y puzzle like it 
is known in (n — d) scattering and may be cured by 3iV-forces. The detailed 
studied 4 He system [10] is unfortunately extremely difficult to describe due 
to the 4 He bound-state and particularly the complex resonance-structure. For 
that reason we investigate the much simpler system (p — 3 He), the 4 Li. In the 
studied energy range up to 30 MeV is a lot of data available. And we have 
restricted ourselves on realistic iViV-interactions without three particle forces 
namely the Argonne-f 14-potential [11] and the Bonn-potential in the form of 
Ref. [12]. 

The paper is organized as follows: In section 2 we give a brief introduction to 
the RRGM. We present in section 3 for each of the interactions (Argonne-f 14 
and Bonn) two model-spaces (a medium and large one) for the used RRGM 
scattering calculation. Then we compare and discuss in section 4 the obtained 
phase-shifts 5, analyzing powers A y and differential cross-section da/dil c . m . 
with experimental data and demonstrate the origin of the /^-problem. 



2 Refined Resonating Group Model with Distortion channels 

2. 1 Refined Resonating Group Model 

This calculation is based on the Refined Resonating Group Model (RRGM) 
[13] with distortion channels [14]. Therefore we briefly summarize the RRGM 
in the following. 

The model employed is restricted to a two fragment description. That means 
the model is not able to describe the three particle breakup reactions. There- 
fore a scattering channel consists of two fragments i — 1,2 with total angular 
momentum Jj which couple to the channel spin S c . The channel spin S c cou- 
ples with the relative orbital angular momentum L TC \ to the total angular 
momentum J of the channel 25c+1 L re ij. The binding energy and the wavefunc- 
tion of one fragment is obtained using the variational principle of Ritz. The 
solution of the Schrodinger-Equation for the scattering problem is won by the 
variational principle of Kohn-Hulthen [15]. 

The Hamilton-operator H for a iV-particle system with a two-body-force is 
given by: 

N N 

H(l,...,N)=Y,Ti+ E Vi, ■ (!) 

1=1 »,J=1 
i<j 

Using momentum conservation, the center-of-mass energy T c m . is separated 
off. The restriction to a two-fragment model allows to formulate a translation- 
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ally invariant Hamilton-operator H' consisting of fragment and relative-motion 
parts. The potential term becomes shortranged by subtracting the Coulomb- 
potential Z x Z 2 e 2 /R with the relative coordinate R between the fragments: 



H'(l, N) = H^l, Ni) + H 2 {N X + 1, N) + T rel + Z X Z 2 e 2 /R 



( 



+ 



£ Vij - Z l Z 2 e 2 /R 



ie{i,...,Ar 1 } 

\j6{JVi+l,.-.,JV} 



(2) 



With the variational principle of Kohn-Hulthen [15] 

5 \H' - E\ *,) - ^a„) = 



(3) 



we determine the solution of the Schrodinger-Equation where the scattering 
wavefunction is denoted by \E^. From the reaction matrix a the scattering 
matrix S is calculated via the Cayley-Transformation: 



S = (]l + ia) (]l-ia) 



-i 



(4) 



The index I in Eq.(3) is the label of the corresponding boundary condition. 
The ansatz for the wavefunction is 



Ik 
rel 



<k=l 



with the antisymmetrization operator 



(5) 



(6) 



where P is the permutation over all particles, n k is the number of channels, 
ty l c is the channel function and ^ei is the relative wavefunction: 

¥ r k el (R)=5 kl F k (R)+a lk G k (R)+J2bi km xUR) • (7) 



Ffc and G k are regular and regularized irregular Coulomb-functions for the 
correct description of the asymptotic behavior of the wavefunction. 5 k i is here 
the Kronecker-Symbol. The \km are square-integrable functions for the de- 
scription of the wavefunction in the interaction region. The variation of the 
a ik and the b\ km yields with Eq.(4) the ^-Matrix: 



S k i = m e 2l5fci 

with phase-shifts 5 k i and channel coupling-strengths rjki- 



(8) 
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The scattering wavefunction (see Eq.(5)) was already explained but without 
the channel part. We look only at one term in Eq.(5). It is convenient to make 
the following ansatz for the channel function \l/ c : 



,(r) 



i J 



(9) 



We use = 1,2) for the translational invariant wavefunction of the ith 

fragment with spin J { . In Eq.(9) square brackets indicate the angular momen- 
tum coupling. In the case of scattering calculation the coordinate r is the 
relative coordinate R. If ^ rc \(r) is a bound-state wavefunction then there will 
be no coupling to the channel spin S c and the coordinate r will be the Jacobi- 
coordinate between the center-of-mass of the both fragments. The fragment 
function cj) J is build of the spin-isospin function S and of the spatial function 

X- 

(10) 



J = E [<# 



LS h^ST 



h,S,a 



The set of inner orbital angular momenta in the fragments is labeled with 
//. The corresponding fragment consists of n nucleons. The spin-isospin func- 
tions S are coupled to good total spin and good isospin. The Clebsch-Gordan- 
coefficients of the coupling of the orbital angular momenta / 1; Z 2 , l n -i to the 
total orbital angular momentum L, the coefficients of the coupling of the spins 
and isospin are coupled to S and T and the coefficients of the super-position 
of the radial dependencies a are expressed by the factor C l J LS . The Xa are the 
square-integrable spatial functions. They consist of nucleon relative functions 



n-1 



k 



k=i 



The nucleon relative function consists from a Gaussian function with width 
parameter (3 a> k and from a Solid-Spherical-Harmonics y [16]: 



Xa,k = eX P {-Pa,k Pi} yi k ,rn k (Pk), 



(12) 



where p k is the Jacobi-coordinate between the (k + l)th nucleon and the 
center-of-mass of the nucleons (1, 2, ...k). 

All spatial functions are in this model finally parametrized with Gaussian 
functions. 



2.2 Distortion Channels 



Distortion channels in the RRGM are unphysical bound-states which can be 
chosen arbitrarily as far as they are linear independent and as long as they 
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have right total angular momentum and total parity. They are allowed to 
violate e.g. the channel spin S c of the scattering channels. Distortion channels 
have the asymptotics of bound-state channels. The idea of including distortion 
channels is to enlarge the model-space in the region of interaction. Therefore 
distortion channels can be used to describe physical structures which can not 
be described by the RRGM or are not included in the scattering channels e.g. 
deformations of the fragments, excited fragments and scattering-states which 
are not yet opened channels. A more detailed description of the approach using 
distortion channels with the RRGM can be found by R. Wolker [14]. 



3 Model-Spaces for (p — 3 He) Scattering 




(a) (b) (c) 

Fig. 1. Panel (a) and (b) show the (p — 3 He)-structures and subfigure (c) shows the 
(d — 2 He)-structure used to describe the (p — 3 He)-scattering-system with L rc \ < 3. 

For the description of the 4 Li-scattering-system with L re \ < 3 and J n = 
0^, 1^, 2^, 3 + we have used the sets of Jacobi-coordinates with orbital an- 
gular momenta as indicated in Fig.l. We have studied this scattering-system 
with the Argonne-wi4- (avu) and the Bonn-potential. For the 3 He-subsystem 
we have used two different model spaces: a medium system (labeled with HeM) 
with 5 basis-vectors and a large 3 He-subsystem (labeled with HeG) with 23 
basis-vectors. The (p — 3 He) -scattering-system using the large 3 He-subsystem 
was studied without F-wave scattering because the F-wave phase-shifts turned 
out very small in the medium system (p — 3 He) -system. Table 1 shows the used 
model-spaces in compact form. The sets of width parameters consisting of n\ 
elements on the first Jacobi-coordinate and ri2 on the second indicated in 
brackets are assigned to spin-isospin-functions with the choice of angular mo- 
menta. The sets of widths are determined by a non-linear optimization [17]. 
In the scattering calculation for the large 3 He-subsystem the smallest width 
on the second coordinate was dropped due to the fact that there was practi- 
cally no improvement of the 3 He wave-function. For the rearrangement channel 
(d — 2 He) we used for the description of the deuteron d three width parameters 
for the S-wave and two for the D-w&ve. The unbound 2 He-system is described 
with an S'-wave using the S'-wave width set of the deuteron. The rearrange- 
ment channel (d — 2 He) opens in the medium system at 5.7 MeV using the 
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Argonne-f i4- and at 6.5 MeV using the Bonn-interaction. For the large model 
these energies are 7.1 MeV (01*14) and 7.6 MeV (Bonn). The binding-energies 
E B , the charge r.m.s. -radii and the D-state probabilities of the models are 
summarized in Tab. 2 and Tab. 3. As expected the D-state probabilities for 
the Argonne-wi4-interaction are appreciably higher than those for the Bonn- 
potential, see table 3. With increasing the model-space the binding-energy 
is improved and also the r.m.s. -radii. The r.m.s. -radii are smaller than the 
experimental ones and they are even smaller for smaller model-spaces. Nat- 
urally one would expect for models which underbind the system to generate 
larger r.m.s. -radii. Here, however, in the smaller model-spaces those width pa- 
rameters, which generate the larger radii, are missing. A basis-vector of the 
4 Li-scattering-system is therefore a spin-isospin function consisting of one neu- 
tron and three protons with a fixed orbital angular momentum configuration 
and the appropriate set of width parameters. The set of distortion channels for 
both models was obtained by using the (d — 2 He*) and (d — 2 He**) channels as 
distortion-channels because they open at energies higher than 30 MeV and are 
therefore out of the studied range of energies. In addition as much as possible 
of the existing basis-vectors have been used to form the maximum of linear 
independent distortion-channels. Every additional distortion channel contains 
therefore only one basis-vector. For the scattering-channel the set of widths 
w rc i was used on the relative coordinate, where for the distortion channels only 
the subset Wdist was taken on that coordinate. 



4 Results and Discussion 

4-1 {p — 3 He) phase- shifts 

In this subsection phase-shifts are diagonal-phase-shifts. Published eigen-phase- 
shifts have been transformed to diagonal-phase-shifts. 

Unfortunately there are no modern theoretical phase-shift calculations at 
higher energies up to -E c . m . = 30 MeV. Only Reichstein et al. [21] and Heiss 
et al. [22] did Resonating Group Model calculations using a much smaller 
model-space and a much simpler AW-interaction in the considered energy 
range. Reichstein et al. used a purely central AW-potential. Heiss et al. used 
a softcore central, spin-orbit and tensor-force potential. Beltramin et al. [23] 
made a separable potential model up to E c m ^ = 10 MeV. We will not com- 
pare their results with ours since their models are much simpler than ours. 
Viviani et al. [24] recently studied the elastic (p — 3 He) -scattering-system us- 
ing the Kohn variational principle with correlated hyperspherical harmonics. 
As interaction Viviani et al. used the av ig-A^/V-interaction [25] and the avig- 
A^V-interaction with the Urbana IX 3A"-potential [26]. Unfortunately they 
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12.95665 


5.134670 


2.947287 


1.342339 


0.821446 


Wdist [fm 


- 2 ] (5) 


0.444741 


0.293900 


0.169016 


0.118524 


0.084300 






0.050011 


0.025737 


0.013852 


0.007143 


0.003852 






0.001857 


0.000973 


0.000562 


0.000277 


0.000101 


Table 1: 


The model-spaces for the 


3 He-system 


is given 


in 



the upper two third of the table. The sets of widths are as- 
signed to the spin-isospin-function and angular momentum 
structures. The number of sets of widths on the first and the 
second Jacobi-coordinates is indicated in brackets. (N is p 
or either n and therefore N ^ N denotes n or p). The sets 
of widths for the relative coordinate in the 4 Li-scattering- 
system for the scattering- w re i and distortion-channels waist 
(underlined) are given in the lower third of the table. 

published only results for very low energies (-E c .m. < 1-69 MeV). Pfitzinger 
et al. [27] did an i?-matrix analysis and an RRGM calculation but only at 
low energies (-E c .m. < 5 MeV). We will not include the results of Pfitzinger 
et al. into the figures. In the following we will compare our results with the 
phase-shift analysis based on experimental data in the desired energy-region of 
Darves-Blanc et al. [28], McSherry et al. [29], Midler et al. [30] and Tombrello 
[31]. 

The phase-shifts obtained by the RRGM calculation using the medium model- 
space are compared with the results of other groups in Fig. 2-4. The S*-wave 
phase-shifts, shown in figure 2, are in very good agreement in the whole range 
of interest with the results obtained by others [24,28,29,30,31]. The S-wave 
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AW-Interaction 


Argonne-fi4 


Bonn 


E B (d) [MeV] 
Eg ax (d) [MeV] 
E^ p (d) [MeV] 


-1.878 
-2.222 


-1.911 
-2.222 [12] 


-2.225 [18] 


£^( 3 HeV 2 +) [MeV] 

^(3 He i/2+) [ MeV j 

£;f ax ( 3 He 1/2+ ) [MeV] 
^ xp ( 3 HeV 2 +) [MeV] 


-5.588 
-6.895 
-6.949 


-6.287 
-7.416 
-7.492 


-7.718 [19] 



Table 2: The binding-energies E-q of the used model-spaces 
for the fragments in the 4 Li-scattering-system are summa- 
rized. For comparison the experimental data (exp) and the 
best result of a RRGM calculation (max) are shown for the 
different studied iViV-interactions Argonne-^u and Bonn. 



iViV-Interaction 


Argonne-vi4 


Bonn 


r r.m. S .(.d) [ fm ] 

r™\(d) [fm] 
rr x m.s.(d) [fm] 


1.80 (5.37%) 
1.98 (6.08%) 


1.77 (4.51%) 
1.98 (4.78%) 


2.13 [20] 


r r M m .,( 3 He 1/2+ ) [H 
r r G m .,( 3 HeV 2 +) [fm] 
« s .( 3 He 1 /2+) [f m ] 
r r ex r g. s .( 3 He 1 /2+) [ fm ] 


1.72 (7.47%) 
1.86 (8.98%) 
1.88 (9.00%) 


1.62 (5.93%) 
1.76 (6.96%) 
1.78 (7.05%) 


1.97 [19] 



Tabic 3: The charge r.m.s. -radii r r . m . s . of the used model- 
spaces for the fragments in the Li-scattering-system are sum- 
marized. For comparison the experimental data (exp) and 
the best result of a RRGM calculation (max) are shown for 
the different studied A?" TV-interactions Argonne-vi4 and Bonn. 
The D-state probability of each model is shown in parenthe- 
ses. 

phase-shifts are all negative due to the underlying Pauli-forbidden states. The 
D- and F-wave phase-shifts are not shown because they are small and they 
are in good agreement with the results of others cited where the .D-wave 
phase-shifts are between —10° and 0° in the whole energy range. The F- 
wave phase-shifts are even smaller than the D-wave phase-shifts by a factor of 
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about 10. For the P-wave phase-shifts, see figure 3, however, the situation is 
different. They are all positive. The RRGM reproduces the over-all behavior 
quite well, but the maximal values are never reached, that means that all P- 
wave phase-shifts are not attractive enough. Furthermore the maxima occur at 
too low energies. The Bonn-interaction reproduces these phase-shifts slightly 
better. Since the Bonn-potential underbinds the 3 He less than the Argonne- 
t>i4-interaction (see Tab. 2). Both forces predict also a too small r.m.s. -radius in 
the studied model-spaces where Argonne-interaction gives a better description 
(see Tab. 3). Generally the Bonn- interaction yields a slightly better description 
of the 3 He subsystem than the Argonne-f i 4 -force and yields therefore to a 
slightly better description (see Tab. 2 and Tab. 3). 

The analysis of McSherry et al. and Tombrello disagree around 5 MeV for the 
1 Pi (Fig. 3 (b)) and 3 Pi (Fig. 3 (c)). From general consideration we consider 
the results of Tombrello more reliable, as they indicate only a weak j-splitting. 
Our results for the 1~ phase-shifts are quite similar to those of Pfitzinger et 
al. [27]. 

But for the 3 P 2 phase-shift analysis (Fig. 3 (d)) the results of the other groups 
are consistent. Here we have a good reproduction at low energies < 3 MeV. 
Then our phase-shift shows clearly not enough attraction. Also Pfitzinger et 
al. [27] reported that they are not able the to reproduce the 3 Po an d the 3 P 2 
phase-shift. Therefore we will concentrate in the following on the 3 P 2 phase- 
shift. 

In Fig. 4 we compare the 3 P 2 phase-shift calculated using the Bonn-interaction 
for the medium and the large model-space with the results of others as indi- 
cated. It turns out that the better 3 He wavefunction does not necessarily lead 
to a better description of the phase-shift although this model reproduces more 
accurately the binding-energy and r.m.s.-radius (see Tab. 2 and Tab. 3). Sim- 
ilar results have been obtained using the Argonne-f 14 -interaction. The overall 
description of the phase-shifts is improved using the large model-space but 
the P-wave phase-shifts show less attraction than in the medium model. In 
the medium model the P-wave phase-shifts drop too fast at higher energies 
because of the too small r.m.s.-radius in comparison to the larger model-space. 



4-2 Analyzing powers A y and differential cross-section dcx/dn c . m . 

To check if the 3 P 2 phase-shift carries the main contribution of the so-called 
Ay-problem we will take now a closer look at analyzing-power A y and the 
differential cross-section da"/dJl c . m .. We compare the analyzing-power A y and 
the differential cross-section at P c .m. = 5.20 MeV and at P c .m. = 22.5 MeV 
with data, namely Birchall et al. [32], McCamis et al. [33], McDonald et al. 
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Fig. 2. The 1 S'o (a) and the 3 5i phase-shifts (b) of (p — 3 He) scattering using 
the medium model-space for the 3 He wave-function for the Argonne-vu (dashed 
line) and the Bonn (solid line) interaction in comparison with experimental and 
theoretical results from other groups is shown: Darves-Blanc et al. [28] (o), McSherry 
et al. [29] (□), Miiller et al. [30] (+), Tombrello [31] (x) and Viviani et al. [24] (o). 

[34] and McSherry et al. [29] for the analyzing-power results and McDonald et 
al. [34] and Harbison et al. [35] for the cross-section results. The error bars are 
increased in comparison to the original ones of McCamis et al. [33] because 
we took the values from the their figures since the data was not tabulated. 

In Fig. 5 we compare the results for the differential cross-section of McDonald 
et al. at £ c . m . = 5.21 MeV and Harbison et al. at P c . m . = 23.0 MeV with 
our calculation in the medium Bonn-model at E cm = 5.20 MeV (solid line) 
and P c .m. — 22.5 MeV (dotted dashed line). The curves in Fig. 5 (dashed line 
and dotted line) correspond to calculations where to the 3 P 2 phase-shift has 
be modified. We have added at E c . m . = 5.20 MeV 10.5° and at E c . m . = 22.5 
MeV 15.0° to reach the experimental 3 P 2 phase-shift in the S'-matrix-element. 
Figure 5 shows that we are now in very good agreement with the experiment 
at P c . m . = 5.20 MeV. But at P c . m . = 22.5 MeV it turns out that a better 
description of the 3 P 2 phase-shift is not enough to explain the experimental 
data. At least we get a good agreement for angles above 110°. 

Since looking on the differential cross-section can only provide a rough overview 
of the quality of the calculation because the differential cross-section is not too 
sensitive to the particle spin. Therefore we study now to the much more sensi- 
tive analyzing-powers (see Fig. 6). We employed the same modifications as for 
the differential cross-section, setting the 3 P 2 phase-shift to its experimental 
value. The calculated analyzing-powers A y demonstrate that 3 P 2 phase-shift 
has the main contribution to the A^-problem. In fact all P-wave phase-shifts 
show too less attraction. Including the 3 P 2 -channel does not the change the 
results. Therefore the tensor-force has no strong contribution to the A y prob- 
lem. 
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Fig. 3. The same as in figure 2 is shown but for the 3 Pq (a), 1 P\ (b), 3 Pi (c) and 
3 P2 (d) phase-shifts. 

4-3 Conclusions 



We have discussed a new RRGM calculation of the 4 Li-scattering system us- 
ing the Argonne-f i4- and the Bonn-interaction and two different model-spaces. 
The calculated phase-shifts have been compared to experimentally and theo- 
retical obtained ones. The overall description of the phase-shifts is in very good 
agreement with the results of others [24,28,29,30,31]. Due to the fact that the 
Bonn-interaction is slightly more attractive than the Argonne-f 14-interaction 
leads to a better reproduction of the experimental data by the Bonn-potential. 
But the P-wave phase-shifts show not enough attraction. This result was also 
obtained by Pfitzinger et al. [27] and Viviani et al. [24]. Pfitzinger et al. and 
Viviani et al. have used the Argonne-f i 8 -interaction with three-particle-forces. 
At low energies the analyzing-power A y can be reproduced by setting the 3 P 2 
phase-shift to the experimental determined one. Pfitzinger et al. [27] pointed 
out that also the 3 Pq phase-shift has to be corrected to get into agreement. 
Unfortunately our model-space is too small to make a statement to this result. 
But at higher energies it turns clearly out that not only the 3 P 2 phase-shift 
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Fig. 4. The 3 i-2 phase-shift is shown using the Bonn interaction for the medium 
(solid line) and the large model (dashed line) for the 3 He wave-function. For com- 
parison results of other experimental and theoretical groups are shown as indicated 
in figure 2. 



has to be corrected to reproduce analyzing-powers and cross-sections. And at 
low energies the contribution of three-particle-forces is negligible since we can 
reproduce all available phase-shifts within our model. The iViV-interactions 
used yield large differences in the P-wave phase-shifts compared to experi- 
mental values. Pfitzinger et al. [27] have demonstrated that new contributions 
to the 3iV-forces acting on P-waves should be considered especially for the 
description of the 3 p2 — 3 Po splitting. Therefore we consider this system well 
suited for further studies of the 3iV-forces. All the realistic iViV-interactions 
studied till now, av' 8 , av u , avi 8 , and Bonn show the same behavior for all the 
P-wave phase shifts. 
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Fig. 5. The differential cross-section dcr/df2 c . m . in comparison with the experiments 
at E c , m , = 5.20 MeV (solid line) and E c , m , = 22.5 MeV (dotted dashed line) is 
shown. The curve where to the 3 i-2 phase-shift 10.5° have be added at £7 c . m . = 5.20 
MeV is the dashed line. And the curve where to the 3 P2 phase-shift 15.0° have been 
added at E c . m . = 22.5 MeV is the dotted line. The experimental data at E c . m . = 5.21 
MeV is taken form McDonald et al. [34] and at E c , m , = 23.0 MeV from Harbison et 
al. [35]. 
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Fig. 6. The analyzing power A y (solid line) is shown in comparison with experi- 
mental results (x) with polarized 3 He at -E c .m. = 5.20 MeV and results of McSherry 
et al. [29] at E c , m , = 5.18 MeV (a), with polarized protons at E c , m , = 5.20 MeV 
and results of McDonald et al. [34] at E CAn . = 5.21 MeV (b), with polarized 3 He at 
E c .m. = 22.5 MeV and results of McCamis et al. [33] (c) and with polarized protons 
at -E c .m. = 22.5 MeV and results of Birchall et al. [32] (d). The curve where to the 
3 P2 phase-shift 10.5° (subfigure (a) and (b)) and 15.0° (subfigure (c) and (d)) have 
been added to reach the experimental results is the dashed line. 
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